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ABSTRACT 

Intranight polarization variability in AGN has not been studied extensively so far. Studying the 
variability in polarization makes it possibly to distinguish between different emission mech- 
anisms. Thus it can help answering the question if intranight variability in radio-loud and 
radio-quiet AGN is of the same or of fundamentally different origin. In this paper we inves- 
tigate intranight polarization variability in AGN. Our sample consists of 28 AGN at low to 
moderate redshifts ( 0.048 < z < 1.036), 12 of which are radio-quiet quasars (RQQs) and 16 
are radio-loud blazars. The subsample of blazars consists of eight flat-spectrum radio-quasars 
(FSRQs) and eight BL Lac objects. Each AGN was observed for a time span of ~ 4 h in the 
R band to measure polarization and variability. Using statistical methods, we determine duty 
cycles for polarized emission and polarization intranight variability. We find clear differences 
between the two samples. A majority of the radio-loud AGN show moderate to high degrees 
of polarization, more than half of them also show variability in polarization. There seems to be 
a dividing line for polarization intranight variability at P ~ 5 per cent over which all objects 
vary in polarization. We did not find clear correlations between the strength of the variability 
and the redshift or degree of polarization. Only two out of 12 radio-quiet quasars show polar- 
ized emission, both at levels of f < 1 per cent. The lack of polarization intranight variability 
in radio-quiet AGN points towards accretion instabilities being the cause for intranight flux 
variability whereas the high duty cycle of polarization variability in radio-loud objects is more 
likely caused by instabilities in the jet or changes of physical conditions in the jet plasma. We 
were able to constrain the timescale of the detected variations to > 4 h. Further studies of 
intranight polarization variability will be necessary to reveal exact physical conditions behind 
this phenomenon. 

Key words: accretion, accretion disks - polarization - BL Lacaertae objects: general - 
quasars : general - galaxies: jets 



1 INTRODUCTION 

AGN are amongst the most powerful sources of radiation in the 
universe. They form giant powerful jets and are thought to be key 
objects in the evolution of galaxies and supermassive black holes. 
Still, we do not yet fully understand the physics behind these ob- 
jects. 

The class of AGN is divided into a wide variety of sub- 
classes. Depending on their luminosity, AGN are divided into the 
low-luminosity Seyfert galaxies and the high-luminosity quasars 
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jVeron-Cettv & Verorj20"o^ . AGN are also classified according to 
their radio-loudness into radio-loud and radio-quiet AGN, where 
the radio-loudness is defined as the ratio between the radio and op- 
tical flux (Kellermann et al. 1989). Additionally, AGN are divided 
into a variety of classe s according to their spe ctral features, polar- 
ization and variability jUrrv & Padovanill 19951) . 

It is an open question what i s the major d iff'eren ce between 
radio-loud and radio-quiet AGN. lUlrich et all ( 1 19971) stated that 
while radio-loud AGN produce powerful collimated jets, such well 
collimated jets are missing in radio-quiet AGN. Differences in the 
spin of the the blac k hole might cause a AGN to be radio-loud or 
radio-quiet (see e.g. lSikora et alJllOOTh . It has also been suggested 
that the viewing angle decides if a quasar is observed to be radio- 
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loud or radio-quiet jKellermann et al.lll989t) . Yet latest studies of 
a radio-quiet quasar ( iBlundell et alj|2003h that showed a relativis- 
tically beamed jet imply that radio-loudness might be a temporal 
phenomenon or that the jets of radio-quiet quasars are simply far 
less powerful. 

In this paper, we will focus on two classes of AGN: radio- 
quiet quasars (RQQs) and blazars. The blazar class is divided into 
two subclasses: BL Lac objects and flat-spectrum radio quasars 
(FSRQs). While in BL Lacs the emission lines are either weak or 
absent, FSRQs show broad emissi on lines. Both typ es are char- 
acterised by a flat radio spectrum ("Ulrich eta lj|l997h. So far, all 
known blazars are radio-loud (Angel & Stock mai] 1980l) . Accord- 



ing to the standard interpretation, blazars are AGN with a j et point- 
ing directly towards the observer jUrrv & Padovanill 19951) . Radia- 
tion from the jet is thought to be highly beamed due to the small 
viewing angle. One of the most striking properties of this class is 
their vi olent variabili ty on all timescal es from hours to decades, 
see e.g. lValtaoia et alj ( l200d) . lstalin et al] ( l2005h . Also, blazars are 
known to show extremely high degrees of polarization and vari- 
ability. In many cases, the most variable objects also sho w the 
highest degrees of polarization ( Angel & Stockman 1980). This 
makes them promising objects for studying intranight polarization 
variability. A de tailed review of polarizati on properties of blazars 
can be found in lAngel & StockmanI 1 19801). B ased on gamma-ray 
data. lGhisellini . Maraschi & TavecchioT 2009h suggested that there 
might be a clear difference between FSRQs and BL Lacs, with FS- 
RQs showing more efficient accretion and thus stronger jets, while 
the accretion power is weak in BL Lacs, resulting in less collimated 
jets and weak line emission. 

In this paper, we will study the variability in polarization 
on time scales of hours, called intranight polarization variability. 
Studying short time scale flux variations in AGN, two different 
terms are used to specify the presence or absence of variability: 
microvariability and intranight variability. While microvariability 
refers to variability on mmag levels, usually on smaller time scales 
dMiller & Noblell996h, intranight variability refers to variability on 
time scales of hours dWagner & Witzelll995f) . Thus, microvariabil- 
ity is defined by the strength of the variations while intranight vari- 
ability is defined by the time scale of the variations. In this study, 
we will use the concept of intranight variability to define the pres- 
ence or absence of variability. 

Intr anight variability have been observed in numerous AGN 
(see e.g. IWagner & Witzel|[T993. iGopal-Krishna et al.ll2000h . The 
pheno menon is mos t prominent in blazars and radio-loud quasars 
(see e.g l Wagner & W itzel 1995). However, even though they show 
less violent variability, also radio-quiet qua sars (RQQs) have 
been found to intra night variability (see e.g. IStalin et all l2005l . 
iGopal-Krishna et al] 2000). Still, the origin of intranight variabil- 
ity is not yet established. 

If one wants to understand the origin of intranight variability 
in AGN, one needs to distinguish between the different mechanisms 
that can cause the variability. A wide variety of theoretical models 
has been proposed to explain variations on small time scales. A way 
to distinguish between these models is using statistical arguments, 
as different mechanisms are fostered by certain conditions. 

Instabilities in the jet are thought to cause violent vari- 
ability in radio-loud AGN, especially in the subclass of blazars 
dWagner & Witzel 1995). Different models exist that explain 
vari ability in jet emission, e.g. t hrough lighthouse effects (see 
e.g. Gopal-Krishna & Wiital 1992h. shock fro nts in jets (see e.g. 
lMarscher&Geaj|l995l. iMarscher et alj|2008h or cer tain physical 
conditions in the jet plasma (see e.g. Bi6mssoi]|l985h . 



The ultraviolet wavelength region in quasars is dominated 
by the rmal emission from the accretion disk, the so-called 'blue 
bump' dGaskelll bOOS'). Thus, variability in quasars at these wave- 
lengths is more likely to be caused by accretion instabilities. This 
dependency can be studied by observing quasars at different red- 
shifts with the same filter and thus at different restframe wave- 
length. However, this makes it possible that any actual redshift de- 
pendency in the variability properties of quasars might be inter- 
preted as wavelength dependency. 

For example, it has been suggested bv lLabita et al] d2008h that 
in the optical bands, RLQs might have a higher contribution from 
the non-thermal component than RQQs. They studied a large sam- 
ple of RQQs and RLQs and were able to model the differences in 
the spectral energy distributions (SED) assuming different spectral 
indexes for the non-thermal component, such that RQQs are dom- 
inated by thermal emission in the optical, while for the RLQs the 
non-thermal component is dominant. Thus, comparing the duty cy- 
cle and strength of INOV in RLQs and RQQs can help answering 
the question if INOV is intrinsic to the non-thermal or thermal com- 
ponent, i.e. the jet or disk emission. 

A phenomenon which has not been studied elaborately so far 
in this context is polarization INOV (PIN OV, referring to var i ability 
in polarization on time-scales of hours). lAndruchow et al] d2005l) 
studied intranight polarization variability in BL Lacs, they found 
very high duty cycles for both polarization and polarization vari- 
ability. Polarized emission in AGN can have different reasons. An 
important origin of polarized emission are relativistic jets. Due to 
the strong magnetic fields in relativistic jets, the emission from the 
jet can be domina ted by synchrotron emission an d is thus highly 
polarized (see e.g. lBiomsson & Blumenthailll982h . Due to unclear 
mechanisms, not only the degree of polarization, but also the po- 
sition angle can change on relatively small time scales (see e.g. 
Holmes et al. 198^. It is not yet clear what causes the changes in 
position angle (see e.g. lBiomssonll 19851) . Also in nearby AGN, po- 
larization due to scattering in the exte nded gas region around the 
nucleus or the disk has been detected dAntonuccill 1 9821) . Different 
models predict more or less l ow polari zed emission from scattering 
in the accretion disk (see e.g . lLaor et: al. 1990, Agol & Bleas 199q|, 
lOoosmann & GaskelllboO? ). However, polarization due to scatter- 
ing is not expected to cause a rotating position angle. Additionally, 
if the distribution of the position angle of polarization is symmetric 
it will cancel out when summing up over the whole object. Thus, 
this kind of polarized emission is unlikely to be observed at high 
redshifts. 

While the observation of flux variability does not give a direct 
hint towards the mechanisms that cause it, polarization variability 
does. Thus, by studying the polarization variability and their link 
to flux variability in different types of AGN, it is possible to find 
out which mechanisms causes the variations. This might be the 
only direct way to examine if intranight variability in radio-loud 
and radio-quiet AGN is caused by the same physical conditions or 
if the origin of intranight variability in radio-loud and radio-quiet 
AGN is fundamentally different. 

Our project combines the usage of polarization measurements 
and samples of radio-loud and radio-quiet AGN at different red- 
shift. We measure the duty cycle of PINOV and it's strength and 
derive statistical differences between the two samples. 

In Section 2 we describe the sample, followed by a description 
of the observations in Section 3. The data reduction is presented in 
Section 4, Section 5 sums up the results of our study, followed by 
Discussion in Section 6 and finally Conclusions in Section 7. 
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Figure 1. Redshift distribution of the quasar sample. RQQs are shown as 
striped surfaces, radio-loud objects (FSRQ and BL Lac type objects) as 
checkered surfaces. 

2 SAMPLE 

We selected the AGN according t o the following criteria from 
the Veron-Cetty & Veron catalogue jVeron- Cettv & Verorill2006h : 
AGN were selected to have a listed apparent magnitude in the R 
band of niR < 16 mag to ensure high enough signal-to-noise with 
exposure times in the range < 60 s at the 2.5-m Nordic Optical Tele- 
scope (NOT). In order to minimise dilution of possible variability 
or polarization due to host galaxy flux, we require a non-detection 
of the host galaxy or Mj^s, > M,„,£.;^„s -♦- 2. For objects where such 
data was not available, we assumed a host galaxy absolute magni- 
tude M* and applied the same criterion. Additionally, we required 
at least one comparison star at similar magnitude within the FOV 
of the ALFOSC calcite plate (circle of ~ 140 "diameter). 

Our sample consists of 28 AGN, 16 of them are radio-loud, 
12 are radio-quiet. Among the radio-loud AGN, eight are BL Lac 
type objects and eight are FSRQs (Flat Spectrum Radio Quasars). 
The redshift range for the sample is 0.048 < z < 1.036 (0.048 < 
z < 1.024 for the radio-quiet objects and 0.138 < z ^ 1.036 for the 
radio-loud objects). Table[T]shows the sample. The redshift distri- 
bution of our sample, separately for the radio-loud and radio-quiet 
AGN is shown in Fig.[T] 



3 OBSERVATIONS 

Observations were carried out at the Nordic Optical Telescope on 
La Palma (Canary Islands, Spain) using ALFOSC, the Andalucia 
Faint Object Spectrograph and Camera, and FAPOL. FAPOL is 
an integral part of the FASU unit used in combination with AL- 
FOSC. The observations were done in visitor mode on 5 separate 
runs in 2006 ~ 2008. For the observations we used ALFOSC in its 
polarimetric mode. The instrument setup consists of a A/2 wave 
plate inside FAPOL and a Calcite in the aperture wheel of AL- 
FOSC that splits the incoming light into an ordinary and an ex- 
traordinary beam. The ALFOSC filter wheel that we used is lo- 
cated after the Calcite in the lightpath, thus the filter does not intro- 
duce errors to the polarization measurements. We used four posi- 
tions of the A/2 wave plate: 0°, 22.5°, 45.0° and 67.5°, correspond- 
ing to a rotation of the plane of the electric field of 0°,45°,90° 
and 135°. We used the fast photometry mode of ALFOSC which 
mak es it possible to quickly rea d out several windows on the chip 
(see lOstensen & Solheimll200oh . Therefore, we only read out the 



two beams for the object and all available comparison stars and 
one additional sky window for each object. As the Calcite used in 
ALFOSC is relatively small, it produces a field-of-view of about 
~ 140"in diameter, thus the number of available comparison stars 
is low. Depending on the object, we were able to observe one 
to three comparison stars. As the single windows are of the size 
13"xl3", we did not use standard sky subtraction for the photom- 
etry using a sky annulus. Instead, we used the sky window that 
was set directly next to the two beams of the corresponding ob- 
jects for sky subtraction. Exposure times were chosen depending 
on the magnitude of the object and the comparison stars, such that 
the read out noise did not dominate the errors for a single mea- 
surement. We used the R band for all observations. The instru- 
mental polarization was checked using zero-polarization standards 
and is negligible. Additionally, possible instrumental polarization 
will be subtracted using the field stars. However, this does not cor- 
rect for possible variations of the instrumental polarization over the 
field-of-view. For other instruments it was found that the instru- 
me ntal polarization is caused by asymmetry in the light path (see 
e.g. lPatat & Romaniellol200^ . As the field-of-view when using the 
Calcite is extremely small, the level of asymmetry in the light path 
is low and therefore, the instrumental polarization should be stable 
within the field-of-view. 



4 DATA REDUCTION 

Data were reduced using ESO-MIDAS Q and rtcn\0. Data in fast 
photometry mode are saved as a series of FITS files containing a 
one dimensional pixel stream and a map of which physical CCD 
pixels are read out. The rtcnv program is used to produce two- 
dimensional FITS images that can be analysed with regular image 
processing software. We also used rtcnv to average the overscan 
pixels in order to determine the bias level, which was subtracted 
from the image. This corrects time variation of the bias level. The 
ALFOSC CCD does not have any spatial structure in the bias level, 
thus subtracting a bias level and not a bias frame is fully sufficient. 

In principle flatfielding is not necessary when using dual beam 
polarimeters as the flatfielding effects cancel out in case the object 
does not move on the chip within a cycle of four exposures that 
form one polarimetric measurement. However, we tested different 
flatfielding methods to check for possible improvements. We dis- 
cuss different flatfielding methods for polarimetry and their advan- 
tages and disadvantages in the following paragraphs. 

We did not use flatfields with the Calcite and retarder due to 
the following reasons: imbalances in ALFOSC cause bending of 
the instrument when moving the telescope. This causes apparent 
movement of the dust specks on the chip. One way to use flatfields 
with Calcite and retarder in the beam without introducing errors 
would be to stop the observations about every 30 min to take dome 
flats and continue observations afterwards. However, this would re- 
quire closing and reopening the dome, which takes about 10 min. 
Additionally, during that time guiding would not be possible. Blind 
tracking might cause the star to move out of the box, causing ad- 
ditional overheads for re-acquisition of the target. Since any gain 
would be lost by the large overheads, this is clearly not a useful op- 
tion. Considering the fact that polarimetry flatfields are not as sim- 
ple as photometry flatfields as they show two overlapping fields, it 
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Table 1. The Q uasar-Sample . Date of Ob servations refers to beginning of night. References for identification: a) IVeron-Cettv & VeronI j2006l) . b) 
IPadovani &Giom mi 1 1995). c) |jackson et"ai] jloO?) 



ID 


Date of 


Observations 




z 




RA 




DEC 


iriR 


type 


Pn 0026+ 1 29 


27 


.09, 


.2006 





.145 


00: 


:29: 


:13, 


.69 


+ 1 3- 16-03 9 




15 


ROO " 


PG 004^+0^9 


26, 


.09, 


.2006 





.384 


00: 


:45: 


:47, 


.27 


+04-10-24 4 

1 . IV/. . ^ 




16 


ROO " 


PROS 10242+0017 


25 


.09, 


.2006 


0, 


.569 


02: 


:42: 


:40, 


.32 


+00-57-27 2 




16 


ROO " 


1 P 05 1 4-0010 


07, 


.02, 


.2007 





.291 


05 


:16: 


:33, 


.50 


-00-77-13 5 




16 


ROO " 


1PS 0647+250 


21 


.01 


.2006 





.203 


06: 


:50: 


:46, 


.60 


+25-03 00 




16 


RI T ac 


S5 0716+714 


20, 


.01 


.2006 


0, 


.300 


07: 


:21 


:53, 


.45 


+71 -20-36 4 




15 


RI T ac 


PKS 07'?5+l 7X 


08, 


.02, 


.2007 


0, 


.424 


07: 


:38: 


:07, 


.39 


+ 17-47 -19 




15 


RI T ac 


QSO B0754+100 


06, 


.02, 


.2007 


0, 


266 


07: 


:57: 


:06, 


.64 


+09:56:34.9 




15 


BL Lac 


IPS 0S06+524 


13, 


.03, 


.2008 


0, 


.138 


08 


:09: 


:49, 


.19 


+52- 1 8-58 4 




15 


RI T ac 


1WGAJ0827. 6+0942 


14, 


.03, 


.2008 


0, 


.260 


08: 


:27: 


:40, 


.10 


+09:42:10.0 




15 


ROO " 


OJ 287 


21 


.01 


.2006 


0, 


.306 


08: 


:54: 


:48, 


.87 


+20 06-30 6 




14 


RI T ac 


OSO B095S+4I4 


23, 


.04, 


.2007 


0, 


.239 


09: 


:56: 


:52, 


.41 


+41:15:22.1 




15 


ROO " 


3C 232 


16, 


.03, 


.2008 





.530 


09: 


:58: 


:20, 


.95 


+32-24-02 3 




16 


FSRQ "'"^ 


PRS 0959+6X5 


17, 


.03, 


.2008 


0, 


.773 


10: 


:03: 


:06, 


.77 


+68-13-16 9 




15 


FSRQ "'^ 


OM 280 


06, 


.02, 


.2007 


0, 


.200 


11 


:50: 


:19, 


.21 


+24:17:53.8 




16 


BL Lac 


OM 295 


15, 


.03, 


.2008 


0, 


.729 


11: 


:59: 


:31 


.83 


+29:14:43.8 




15 


FSRQ" 


ON 325 


13, 


.03, 


.2008 


0, 


.237 


12: 


:17: 


:52, 


.08 


+30:07:00.6 




14 


FSRQ "-^ 


3C273 


08, 


.02, 


.2007 


0, 


.158 


12: 


:29: 


:06, 


.70 


+02:03:08.6 




12 


FSRQ "-^ 


PG 1246+586 


18, 


.03, 


.2008 


0, 


.847 


12: 


:48: 


:18, 


.78 


+58:20:28.7 




16 


BL Lac 


3C 279 


24, 


.04, 


.2007 


0, 


.536 


12: 


:56: 


:11, 


.17 


-05:47:21.5 




12 


FSRQ" 


PG 1254+047 


25, 


.04, 


.2007 


1 


.024 


12: 


:56: 


:59, 


.93 


+04:27:34.4 




16 


RQQ" 


PG 1259+593 


14 


.03, 


.2008 


0, 


.472 


13: 


:01 


:12, 


.93 


+59:02:06.7 




15 


RQQ" 


PG 1307+085 


07, 


.02, 


.2007 


0, 


.155 


13: 


:09: 


:47, 


.03 


+08:19:49.3 




16 


RQQ" 


CSO 873 


18, 


.03, 


.2008 


1 


.014 


13: 


:19: 


:56, 


.23 


+27:28:08.2 




15 


RQQ" 


PKS 1656+053 


25, 


.04, 


.2007 


0, 


.879 


16: 


:58: 


:33, 


.45 


+05:15:16.4 




15 


FSRQ 


PG 1700+518 


24, 


.04, 


.2007 


0, 


.292 


17: 


:01 


:24, 


.87 


+51:49:21.0 




14 


RQQ" 


PG 2112+059 


25, 


.09, 


.2006 


0, 


.457 


21: 


:14 


:52, 


.57 


+06:07:42.5 




16 


RQQ" 


CTA 102 


27, 


.09, 


.2006 


1, 


.036 


22: 


:32: 


:36, 


.41 


+ 11:43:50.9 




16 


FSRQ "-^ 



should additionally be kept in mind that these corrections are non- 
trivial. 

On the other hand, it is unclear how well flatfields with only 
the filter in the beam work. It is unclear how well they correct for 
large-scale variations as the lightpath changes significantly when 
inserting the calcite and retarder. We performed tests concerning 
this matter and found that this kind of flatfields correct at least for 
the very large scale variations, i.e. illumination patterns. However, 
they do definitely not correct for dust specks due to the strong dif- 
ferences in the lightpath. On the other hand, pixel to pixel varia- 
tions, i.e. differences in the gain between the pixels, should be cor- 
rected by these kind of flatfields. We performed tests to check this 
assumption. We checked the scatter of the measured normalised 
Stokes parameters for flatfielded and non-flatfielded data. Results 
where inconclusive, some targets showed slight increase, some 
slight decrease in scatter when using flatfielded data. It seems that 
either the gain differences between the pixels or the exact position- 
ing of the windows on the chip are not completely stable. It is also 
possible that the gain differences are sensitive to the polarization of 
the detected light. Thus, we decided not to flatfield the data. 

As we do not flatfield our data, we expect problems with the 
photometry. Even we use guiding for our observations, the objects 
move on the chip during the 4 hours of the observations, usually on 
scales of 1-5 pixels. This changes the influences of dust specks and 
thus can cause false variability, especially if we consider that mi- 
crovariability might occur on scales of mmag. We performed tests 
on the two objects that have three comparison stars (IE 0514-0030 
and CTA 102). In both cases variability is detected in all magnitude 
measurements. Thus, we conclude that photometry from polarimet- 



ric data is not reliable if proper flatfielding cannot be performed. 
This can be a serious problem for polarimetry as well. In some 
cases, not the degree of polarization but the polarized flux might 
be of interest a. If reliable photometry is not available, errors in 
polarized flux occur. 

The flux was measured using aperture photometry. The flux 
in each beam was measured using circular apertures, the sky was 
subtracted using the corresponding sky region. We skipped the first 
and last ten pixels of the sky regions due to read-out problems at 
the beginning of each window showing as an increasing bias level. 
We tested two methods for the sky determination: the median value 
and the mode value determined as follows: 

mode = 3 X median - 2 X mean (1) 

The mode value should avoid wrong sky determination in 
presence of faint objects in the sky field. These might not be visible 
when inspecting the single frames, but might show in a skewed dis- 
tribution of sky values. We ran the data reduction with the two dif- 
ferent sky value determination methods and checked for differences 
in the scatter of the normalised Stokes parameters. The mode sky 
determination method yielded better results. We measured a clear 
decrease in the scatter for more than 95 per cent of the normalised 
Stokes parameters. We thus use the mode sky determination for all 
our data. 

The best aperture was determined by running the reduction for 
a wide range of apertures, we then choose the aperture that yields 
the lowest scatter in the normalised Stokes parameters and P,, for 
the main target. In case of bad data, i.e. clouds passing, almost fully 
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diluting the object, we reject these bad data by setting a hmit for 
the signal-to-noise-ratio. This method yields moderately to signifi- 
cantly bigger apertures than optimisation using the signal-to-noise- 
ratio in the target. As our study is clearly focused on polarimetry 
and not photometry we decided to optimise all data reduction for 
this purpose. 

The degree of polarization and position angle of polarization 
were calculated as follows: 

Qi = — (2) 

Qm = Go- + 222.5" + 245" + 267.5" 0) 
^ 2o° - 245" 

Qm 

„ 222.5" - 267.5" 



P = ( ^Pl + Pjj X 100% (6) 

PA = 28.7 X atanl (p,,, P,) (7) 

With Ei and O, being the flux in the extraordinary and ordinary 
beam for the angles ; = 0°, 22.5°, 45°, 67.5°, respectively. For the 
error calculations, we included shot and read-out noise. Errors in 
polarization were calculated using error propagation for the above 
formulae. 

We subtracted background polarization using the field stars. 
This assumes that both the quasars and the corresponding field stars 
are affected by the same background polarization. This assumption 
was tested by comparing the polarization of different field stars 
within one field, if available. In most of the cases all field stars 
within one field yielded the same degree of polarization and posi- 
tion angle within error limits. However, this estimate of background 
polarization yields only a lower limit. Additional background po- 
larization caused by dust aligned with the Galactic magnetic field 
might be present. This might for example be the case if the com- 
parison stars are nearby objects and are affected by less absorption 
and thus less background polarization. We will comment on this 
in the result section for the objects for which the determination of 
the background polarization was critical. Background polarization 
was subtraced vectorially, i.e. in f and Py with P^./y being the nor- 
malised Stokes parameters. 

Data averaging for the polarization and position angle was 
performed by averaging P^ and Py and calculating the averaged 
polarization from the background subtracted values. Errors for av- 
eraged polarization were calculated from the scatter of P, and 
Py using error propagation. Additionally, we considered the bias 
due to non-gaussian distrib ution of the degree of polarization from 
ISimmons & Steward ( Il985h . The unbiased degree of polarization 
was calculated as follows: 

We consider all objects for which the unbiased degree of polariza- 
tion P unbiased IS cousistcut with zero polarization to be unpolarized. 
For degrees of polarization P < 2 per cent we give an estimate of 
the asymmetric behaviour of the error distributions derived from 
ISimmons & Steward 11983) . 



We apply zero polarization angle correction using measure- 
ments of high polarization standards. This correction was only ap- 
plied to the final averaged data to allow comparison with other data. 
The calibrated position angle of polarization PA is listed in Table 

m 

Variability detection is challenging. Normal -tests compare 
the scatter with the mean calculated error. This poses the problem 
that faulty error calculations cause faulty detections of variability. 
For the photometry we use error estimates that are standard for 
aperture photometry as used in IRAlB however other sources of 
error than shot and read-out noise are expected. For pol arimetry 
the case is even more severe. IPatat & Romaniellol ( l2006h studied 
elaborately the influences of systematic errors in polarimetric mea- 
surements. These systematic errors are found to depend on the in- 
strument, the telescope and the number of retarder plate positions 
used. Other causes of systematic errors may exist. Additionally, the 
error propagation used only considers shot and read-out noise. In- 
complete error estimates for the photometry thus also cause under- 
estimation of the errors for polarimetry. As the detection of vari- 
ability is the main goal of this study, normal -tests as described 
above are not sufficient. 

We use the analysis of variance (ANOVA) to detect PINOV. 
This statistical test is robust and does not rely on error estimations 
but uses the scatter of the data to detect variabilit y. Additionally, 
ANOVA has been used by other authors (see e.g. Ide Diego et al.l 
ll998llStalinetal.ll2005h that studied INOV in RLQs and RQQs. 

We perform a one way ANOVA using the stats package from 
SciPjQ. We divide the data into bins of ten data points. In our re- 
sult table we give the probability that the data in the different bins 
are not drawn from the same parent population, i.e. the probability 
Pyar that the object is variable. For p,,,,- > 95 per cent in either P, 
or Py the object is considered as variable. We perform ANOVA on 
the values of Pj^ and Py as the degree of polarization shows non- 
gaussian distribution at low degrees of polarization. ANOVA ex- 
plicitly assumes Gaussian distributions, thus it cannot be used for 
the degree of polarization. 

Additionally, we inspect all data by eye. Some of the data 
where ANOVA clearly favours variability show weird variability 
patterns, i.e. after averaging and plotting the data, we see a single 
heavy outlier and otherwise stable normalised Stokes parameters. 
For such cases we perform the following reliability check. We re- 
move the datapoints that cause the single heavy outlier and rerun 
the ANOVA on the rest of the data. If the second ANOVA does not 
detect variability, the variability is 'dubious'. 

We present the results of the statistical tests and other impor- 
tant properties, i.e. the average polarization, the number of com- 
parison stars and the exposure time in Table |2] The plots of the 
polarization variability are shown in the Fig.l2lto[m 



5 RESULTS 

In the following paragraphs we discuss the results for selected tar- 
gets. We do not present the results for targets that are unpolarized 
and do not show variable polarization. Results for all targets can be 
found in Table|2] Plots of normalised Stokes parameters and Py 
for all targets are presented in Figures[2ltoll II 

lES 0647-1-250 is a BL Lac type object at redshift z = 0.203. 

^ Image Reduction and Analysis Facility, distributed by NOAO, operated 
by the AURA, Inc. under agreement with the NSF 
"* see ,http://www.scipy.org] 
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Figure 2. Polarization intranight variability. We applied a signal to noise 
limit for each object. The resulting data were binned with a bin size of ten. 
Data were averaged. We plot the mean and standard enor. Dashed horizon- 
tal lined indicate /"j/, = 0. We plotted normalised Stokes parameters PyjPy 
against time in minutes. Background polarization is not subtracted. Plots 
continue in FiglSlto fTT] 



Figure 3. Polarization intranight variability, continued. For caption see Fig. 

HI 



During our observations this object showed polarization of P ~ 5 
per cent after subtraction of background polarization. However, 
within the four hours of observations our statistical tests as well as 
eyeball inspection of the data do not show variations in the degree 
of polarization. The object is considered non-variable in polariza- 
tion. 

S5 0716+714 is a BL Lac type object at redshift z ~ 0.3 
( iNilsson et al.l2008h . We detect polarization of P ~ 6 per cent after 
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Figure 4. Polarization intranight variability, continued. For caption see Fig. 



Figure 5. Polarization intraniglit variability, continued. For caption see Fig. 

HI 



subtraction of background polarization. One of the Stokes parame- 
ters {Py) shows variability in form of a clear decrease of the degree 
of polarization. The other Stokes parameter (P,) shows variability, 
however, on a much lower level. 

PKS 0735+178 is a BL Lac type object at redshift z = 0.424. 
This object is polarized (P ~ 4 per cent) after background subtrac- 
tion. Statistical tests show variability in Py but not in P, . 

QSO B0754-I-100 is a BL Lac type object at redshift z = 



0.266. The object is highly polarized (P ~ 1 1 per cent) after back- 
ground subtraction. Statistical tests imply variability in both nor- 
malised Stokes parameters. The data show a slow increase in Py, 
accompanied by a simultaneous decrease in P, . We consider the 
object variable in polarization. 

lES 0806-H524 is a BL Lac type object at redshift z = 0.138. 
After subtraction of background polarization, we detect polarized 
emission (P ~ 2 per cent). Statistical tests imply violent variability 
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Figure 6. Polarization intranight variability, continued. For caption see Fig. 



Figure 7. Polarization intranight variability, continued. For caption see Fig. 

HI 



in Fv- From eyeball inspection the variability pattern seems to show 
violent variability on very short timescales. 

1WGAJ0827.6-I-0942 is a RQQ at redshift z = 0.260. For 
this object, one of the two comparison stars showed high levels 
of polarization (~10%), we thus rejected this star for subtraction 
of background polarization. We used the other comparison star for 
background subtraction, this yields bigger errors than using two 
comparison stars. However, we can safely assume that a degree of 



polarization as high as measured in the first comparison star can not 
be due to background polarization, especially as the second com- 
parison star shows low levels of polarization. After subtraction of 
background polarization, we do not detect polarized emission. The 
object does not show variability in polarization. 

OJ 287 is a BL Lac type object at redshift z = 0.306. With 
a degree of polarization of P ~ 16 per cent, this is one the 
most highly polarized objects in our sample. Statistical tests im- 
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Figure 8. Polaiization intranight variability, continued. For caption see Fig. 
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ply variability in both normalised Stokes parameters. This result 
has been verified by eyeball inspection of the data. We see an 
decrease in Py, accompanied by an increase in P,-. OJ287 is one 
of the most well observed blazars. The polarization and variabil- 
ity beha viour has been extensively studied by numerous groups 



ity beha viour has been extensively studied by numerous groups 
(see e.g. iHolmes et alj|l98 1 ISillanpadll"99ll iTakalo et al.N 19921 
iHagen-Thoml ll980llValtaoiaet alj|200a) . During our observations 



Figure 9. Polarization intranight variability, continued. For caption see Fig. 

HI 



in Januar y 2006, this object w as in the decrease phase of a major 
outburst Jvaltonen et"al]|2008l) . 

FBS 0959+685 is a FSRQ at redshift z = 0.773. After sub- 
tracting background polarization, we detect polarized emission 
(P ~ 2.4 per cent). However, the object does not show variability 
in polarization. 

OM 280 is a BL Lac type object at redshift z = 0.200. After 
subtraction of background polarization, we detect minimal polar- 



10 C. Villforth et al. 



PKS 1656+053 



PG 2112+059 



-0.4 - 

-o.eL 



h<i|ilhi'ii/'i 












1 


1 


,1 





100 150 

Time [min] 



-0.2 - 

„ -0.4 - 

^ -0.6 - 

^ -0.8- 

-1.0 - 

-1.2 — 
0,40 — 

0.35 - 
0.30 - 
" 0.25- 
^ 0.20 - 
0.15- 
0.10 - 
0.05 - 

o.ooL 



i * \ 



50 100 150 200 250 

Time [min] 



PG 1700+518 



i:,""i"r'iV'ii 

-0.5 ' ' ' ^ 

g 0.4 - I 



* \ , 



. ' I 



II' I 



100 150 200 

Time [min] 



100 150 200 

Time [min] 



Figure 10. Polarization intraniglit variability, continued. For caption see 

Fig.m 



Figure 11. Polarization intranight variability, continued. For caption see 
Fig.lH 



ization of P = 1.25^15 '^ per cent. No variability in polarization is 
detected. 

OM 295 is a FSRQ at redshift z = 0.729. This object shows 
one of the highest degrees of polarization in our sample (P ~ 15 per 
cent after background polarization subtraction). Violent variability 
is detected in both Stokes parameters, both show a clear decrease 
in degree of polarization. 

ON 325 is a FSRQ at redshift z = 0.237. After subtracting 
background polarization we measure a degree of polarization of 
P ~ 10 per cent. Statistical methods imply clear variability in both 
Stokes parameters. Eyeball inspection show that P,. shows a very 
clear increase over the observing time, whereas the variability in Py 
shows only one 'bump'after about ISOmin of our data stream. 

3C 273 is a FSRQ at redshift z = 0.158. After subtracting 
background polarization, we do not detect polarized emission in 
this object. Statistical tests imply variability in both Stokes param- 
eters. Eyeball inspection verify that both P,. and P,, indeed seem 
to show variability. This might imply that the objects showed very 
low variable polarization that cancelled out after averaging over 
the whole time span of observations. We also detected variability 
in the polarization of both comparison stars. However, the variabil- 
ity in the comparison star was only caused by the first three data 
points, whereas the polarization variability in the object also shows 
when the first three data bins are removed. This seems to exclude 
variability in background polarization. 



PG 1246+586 is a BL Lac type object at z = 0.847. We detect 
polarized emission of P ~ 5 per cent after background polarization. 
No variability in polarization is detected. 

3C 279 is a FSRQ at z = 0.536. In our sample, this is the object 
that shows the highest degree of polarization with a value of P ~ 26 
per cent after background polarization subtraction. Statistical tests 
imply variability in Pj, showing in an increase, whereas P,. does 
not show variability. 

PG 1254+047 is a RQQ at z = 1.024. This object was ob- 
served bv'Gopal -Krishna et al.l l ll993h . lGopal-Krishna et al] feOOOl) 
and [Jang (2005) and none of those detected microvariability. Af- 
ter subtracting background polarization, we detect low-level polar- 
ized emission of P = 0.85^0 3^ per cent. Please note that this might 
well be background polarization, in case the comparison stars are 
at close vicinity inside the Milky Way and the object if aff'ected by 
additional background polarization. No variability in polarization 
is detected. 

PG 1259+593 is a RQQ at redshift z = 0.472. After subtrac- 
tion of background polarization, the measured polarization is con- 
sistent with zero polarization. Statistical test imply variability in P,.. 
Indeed, the polarization seems to oscillate. Variable polarization is 
not detected in the comparison stars. 

PG 130 7+085 is a RQQ at re dshift z = 0.155. The object was 
observed by Ijang & Milled ( 1 19951) . they did not detect microvari- 
ability. We do not detect polarized emission. The variability detec- 
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tion in is most likely caused by a single outlier that is visible in 
the plots. As a reliability check, we remove the ten datapoints pro- 
ducing the outlier and rerun the ANOVA on the rest of the data. In 
this second ANOVA test the probability for variability dropped to 
Pvar = 89.6%. Thus the variability is caused by a single point and 
hence dubious. 

PG 170 0+518 is a RQQ at redshift z = 0.292. This object was 
observed bv lCarini et ai]l l2007l) , no microvariability was detected. 
After subtracting background polarization, we measure low-level 
polarization (P = 0.87^55^.J per cent). However, in case of compari- 
son stars in close vicinity inside the Milky Way, this might well be 
background polarization. No variability in polarization is detected. 

CTA 102 is a FSRQ at z = 1.036. After subtracing back- 
ground polarization, we measure polarized emission of P ~ 3 per 
cent in this object. Statistical tests imply variability in both and 
Py, showing a 'bumpy decrease'of polarization over time in both 
values. 

Overall, we observed 28 objects, 12 RQQs, eight FSRQs and 
eight BL Lacs. We detected polarized emission in 15/28 objects. 
The degree of polarization ranged from as low as P = 0.85^q3* 
per cent to as high as P ~ 26 per cent. Two RQQs (PG 1254+047, 
PG1700+518) showed polarized emission. For the radio-loud ob- 
jects, five FSRQs and eight BL Lacs showed polarized emission. 
We detected variability in polarization in eleven objects. In two 
cases, we detect variable polarization in unpolarized objects (3C 
273, PG1259+593). Nine objects are polarized and variable in po- 
larization. All of these objects are radio-loud, four of them being 
FSRQs, five of them being BL Lacs. Most of the variability shows 
in a increase or decrease of P^iy over the time span of our obser- 
vations. In some cases, this trend is modulated by variations on 
shorter timescales. We will discuss the results further in the follow- 
ing section. See Table|3]for a summary of the sample statistics. 



6 DISCUSSION 

In this section we will sum up and discuss our results, we will calcu- 
late duty cycles for PINOV and polarized emission for our samples. 
The duty cycle describes the percentage of time a system is in a 
certain state, i.e. in our case the time that an AGN shows polarized 
emission or PINOV. It should be noted though that we observed 
our objects for only 4 h and not a full night. Thus, we might have 
missed variability or polarized emission in some objects as we did 
not observe them long enough. We will systematically compare the 
samples of radio-loud and radio-quiet AGN. We will also compare 
our results to those of other authors and discuss possible theoretical 
explanations for the observed behaviour. 



6.1 Statistical properties and violent variability 

It indeed seems clear from our study that radio-loud and radio- 
quiet AGN show distinctly different behaviour regarding intranight 
polarization variability. 13 of the 16 observed radio-loud objects 
showed polarized emission, ranging from as low as P = 1 .25 ± 0.85 
per cent to as high as P ~ 26 per cent. The three radio-loud AGN 
that did not exhibit polarized emission were FSRQs. More than half 
of the radio-loud AGN show PINOV. 

Only two of the observed 12 RQQs showed polarized emis- 
sion. In both cases the measured polarization was extremely low 
and the detection might represent residual background polariza- 
tion. Only one RQQ, PG 1259+593, shows variable polarization, 
however this objects is classified as unpolarized. Such a behaviour 
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Figure 12. Degree of polarization plotted against redshift for different sam- 
ples. 



might either be explained by a flickering low level polarization or 
might be 'false'polarization variability. For example, if there is un- 
derlying constant but low polarized emission from e.g. scattering in 
the host galaxy or the dust torus, non-polarized varying flux from 
another component (e.g. the accretion disk) can cause variability in 

The BL Lac type objects build the most extreme sample, all 
of them show polarized emission from P ~ 1.3 per cent to as much 
as P ~ 16 per cent. Five of the eight observed BL Lac type objects 
showed variability within an observing span of about ~ 4 h. Three 
of the eight observed BL Lac type objects did not show variabil- 
ity in polarization. All objects with degrees of polarization greater 
than 5 per cent are variable in polarization. However, also some ob- 
jects showing degrees of polarization lower than 5 per cents show 
variability in polarization. 

The FSRQs conform a less homogeneous sample. Three of 
the observed eight FSRQs did not show polarized emission during 
the observations. On the other hand the object with the highest de- 
gree of polarization detected in this study, 3C 279 with a degree 
of polarization as high as ~ 26 per cent is a FSRQ. With degrees 
of polarization from P ~ 2.5 per cent to P ~ 26 per cent these 
values span a wider range than for the BL Lac type objects. We ob- 
served one FSRQ that was unpolarized but variable in polarization 
(3C 273). This can be explained in the same manner as mentioned 
above for the RQQ showing the same behaviour. 

However, the differences in the range of measured polarization 
between the BL Lacs and the FSRQs might well be coincidental. As 
the polarization is variable, we might have observed the objects in 
a very different state at another date. 

The observed variability appears in mostly slowly increasing 
or decreasing P, or P,,, sometimes modulated with faster vari- 
ability. Thus it seems that we clearly undersample the variation 
timescales, we can only give a lower limit for the timescales. 

A plot of polarization variability and polarization properties 
against redshift can be seen in Fig. [12] The plot shows that all ob- 
jects showing degree of polarization higher than 5 per cent show 
variability in polarization. However, we detect objects with degrees 
of polarization lower than 5 per cent that show polarization vari- 
ability. We conclude that there might be a correlation between the 
degree of polarization and PINOV in the sense that all objects with 
P > 5 per cent showing variability, but not all variable objects 
showing degrees of polarization > 5 per cent. 
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Table 2. Result table. Description of columns: ID: Object name; z: redshift; #: number of comparison stars observed; exp. time/cycle time: exposure time for 
one angle, cycle time is exposure time per angle plus read-out time; best aperture: radius of best aperture for data reduction in arc seconds; pvar [%]: probability 
for variability in P^^ and Py separately, meaning of flags: N=not variable, Y= variable, D=dubious, see text for discussion; V/RV: redshift corrected violence 
and relative violence of variability, see text for explanation; P: degree of polarization in %; PA corrected position angle of degree of polarization; type: object 
type. 
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56.25 /N 


4.79 /N 




4.73±0.31 


123 


BLLac 


3C 279 


0.536 


1 


21.7 (25.0) 


1.90 


98.92 /Y 


13.54 /N 


1.23/0.05 


25.85±0.35 


93 


FSRQ 


PG 1254+047 


1.024 


2 


15.6 (20.0) 


1.90 


66.16/N 


36.92 /N 




0.85:«:^6 


170 


RQQ 


PG 1259+593 


0.472 


2 


50.6 (55.0) 


2.28 


95.24 /Y 


72.12 /N 






RQQ 


PG 1307+085 


0.155 


2 


25.4 (30.0) 


1.90 


>99/D 


22.87 / N 








RQQ 


CSO 873 


1.014 


1 


26.8 (30.0) 


1.52 


9.93 / N 


66.97 / N 








RQQ 


PKS 1656+053 


0.879 


2 


15.6 (20.0) 


1.52 


48.91 /N 


88.24 / N 








FSRQ 


PG 1700+518 


0.292 


1 


11.6(15.0) 


1.71 


81.54 /N 


74.42 / N 




0.87!«:^l 


48 


RQQ 


PG 2112+059 


0.457 


2 


56.1 (60.0) 


2.09 


18.33 /N 


92.12 /N 








RQQ 


OTA 102 


1.036 


3 


25.0 (30.0) 


1.33 


96.73 / Y 


>99/Y 


4.81/1.87 


2.57±0.15 


6 


FSRQ 



Table 3. Duty Cycles of polarized emission and PINOV for different samples. Percentages refer to the fraction of sources 
being in the respective state. 



Sample 




P:yes 


P:no 


P:yes,var 


P:yes,not var 


P:no,var 


all 


whole sample 


(#) 


15 


13 


9 


6 


2 


28 


whole sample 


(%) 


54 


46 


32 


21 


7 




RL(all) 


(#) 


13 


3 


9 


4 


1 


16 


RL(all) 


(%) 


81 


19 


56 


25 


6 




RL (FSRQ) 


(#) 


5 


3 


4 


1 


1 


8 


RL (FSRQ) 


(%) 


63 


38 


50 


13 


13 




RL (BL Lac) 


(#) 


8 





5 


3 





8 


RL (BL Lac) 


(%) 


100 





62 


38 







RQQ 


(#) 


2 


10 





2 


1 


12 


RQQ 


(%) 


17 


83 





17 


8 





To check the probable dependency of variability on degree of 
polarization, we study an estimate of the intensity of the variations: 
the violence V and the relative violence RV. We define the violence 
of the polarization variations as follows. First, we calculate the 
standard deviation of the binned datasct (bin size 10) for and Py, 
CTioiai- Additionally, we estimate the measurement error as the mean 
standard error of the binned data (bin size 10), 5 . In principle this 



kind of calculation could also be done with the standard deviation 
of the unbinned data as (Ttotai and the mean error of the unbinned 
data as 5 . Both possibilities have their advantages and disadvan- 
tages. Considering a;„,ai, this value is much better constrained for 
the unbinned data as the number of datapoints is ten times the num- 
ber of datapoints for the binned data. On the other hand, S is well 
defined for the binned data as it represents an empirical estimate of 
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the error, whereas S for the unbinned data is affected by the fact 
that we only consider shot and read-out noise and other sources 
of error are not considered. We already discussed in Section 4 that 
systematic errors are i mmanent to dual beam pola rimeters as used 
for this study (see e.g. IPatat & Romaniellcll2003) . We argue that 
the presence of systematic errors poses a more severe problem than 
the possibly not well constrained standard deviation when using 
binned data. Thus, we decide to use the binned data to estimate the 
violence V and relative violence RV . 

From the two values calculated as described above, we esti- 
mate the real variability range tr in P^iy as: 



Redshift dependence of violence V 



• CBL Lac 

T TFSRQ 



(9) 



From these values we estimate the range of the variation r, repre- 
senting the FWHM of the distribution of the degree of polarization 
P in polarization corrected for measurement errors as follows; 



r = 2.35 X ^(r(P,y- + o-{Py)^ 
Inserting the expressions for a{P^iy) this yields: 



(10) 



Dependence of violence V on degree of polarization P 



• CBL Lac 

T TFSRQ 



r = 2.35 X ^cr,„aP.f- -S{P,y + cr^UPy)- - 



(11) 



Using this value and the timespan of observations Af and the red- 
shift z we calculate the redshift corrected violence V: 



V 



At 



x(l +z) 



(12) 



In the case of lES 0806-1-524, cr,„,„i < S and thus the violence 
V is not defined. We interpret this as a sign that lES 0806+524 
is marginally variable. As discussed above, the standard deviation 
Ciomi is not well constrained for binned data. Thus it is possible that 

C total < S . 

We plot the violence V against both the redshift z and the de- 
gree of polarization P in Fig. [T3] We do not find any dependency 
between the degree of polarization and the violence. This somewhat 
stands in opposition to the finding that all objects with a degree of 
polarization P > 5 per cent show variability. One would expect that 
if variability is linked to high degrees of polarization, there would 
be a dependency between the intensity of variability and the degree 
of polarization. As we do not find such a dependency, this may ei- 
ther mean that low number statistics conceal possible dependencies 
or that the finding itself is caused by low number statistics. On the 
other hand it may also mean that there is a lower limit in degree of 
polarization and variability for a certain phenomenon. However, it 
is unclear what might cause such a dividing line. A possible cause 
for the seeming dividing line for PINOV might be the fact that de- 
tecting variability in low polarized objects is much more challeng- 
ing. If we consider a highly polarized object with P = 25 per cent a 
comparatively change proportional to the degree of polarization is 
easily detectable, whereas the same change proportional to the de- 
gree of polarization will not be detectable with the same error range 
at « 1 per cent. Bigger samples will be necessary to answer this 
question. 

We do not see a correlation between the violence V and the 
redshift z with the violence increasing with redshift. 

However, the violence might not be the best possible measure 
for the intensity of the variability. Thus we additionally define the 
relative violence RV as follows: 



Figure 13. Dependency of Violence V on both redshift (Upper Panel) and 
degree of polarization (Lower Panel). Violence is in units of lO^^min"'. 
Non- variable polarized sources are also shown in the plot, for those we set 
the violence to V = 0. However, we do not plot the values for the two 
possibly RQQs. 



RV 



{P) X Af 



x(l +z) 



(13) 



This value might be a better measure as it is normalised with the 
mean degree of polarization {P), calculated from the mean of the 
normalised Stokes parameters Pxiv- Whilst the violence would 
yield the same value for a change from P = 25-26 per cent and 
a change from P = 1-2 per cent, the relative violence considers 
the comparative change in degree of polarization. However, calcu- 
lating the relative violence RV in this way raises the problem that 
the correlation between RV and P is somewhat implicit as RV is 
calculated using P. We still show the same plots as for the violence 
Vin Fig. [14] 

As expected we see a tendency in the correlation between the 
relative violence RV and the degree of polarization, with higher de- 
grees of polarization showing lower relative violence. As discussed, 
this tendency might be implicit as RV is calculated using P. 

There seems to be a weak positive correlation between the rel- 
ative violence RV and the redshift. Still, with only nine datapoints 
we are clearly limited by small number statistics. A possible depen- 
dency on redshift might be caused by three different conditions. On 
the one hand there might be a real redshift dependency of variabil- 
ity intensity, on the other hand, the observed effect might be caused 
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Table 4. Spearman rank-order correlation coefficient 
p and corresponding p-value (in brackets) for differ- 
ent combinations of the variables violence V, relative 
violence RV, degree of polarization P and redshift z- 



P 



V -0.17 (0.69) 0.12(0.78) 
RV -0.67(0.07) 0.17(0.69) 



Redshift dependence of relative vioience RV 



• CBL Lac 

TTFSRQ 



▼ • 



Dependence of relative vioience RV on degree of polarization P 



• CBL Lac 

T TFSRQ 



0.0 - •* • 



the corresponding p-value using the statistic package of SciPj[f[ p 
denotes the probability that the measured correlation coefficient or 
a higher correlation coefficient originates from random data. How- 
ever, it should be noted that p can not be calculated reliably for 
small datasets with <500 data points, we still calculate p, but the 
reader should keep in my mind that the values merely represents 
estimates. 

For the relations between the variables V and z, V and P, RV 
and z, the Spearman rank-order test indicates small correlations. 
Correlations of both violence parameters and the redshift are posi- 
tive, whereas the correlations between the violence parameters and 
the degree of polarization are negative. However, as stated above, 
we are clearly limited by small number statistics. This fact also 
shows in the p-values, for the relations z-V, z-RV and P~RV the 
probability that the measured weak correlation is caused by random 
data is > 50 per cent. Thus we can safely say that we did not detect 
any correlations for these values. For the parameters RV and P, the 
Spearman test indicates a strong negative correlation, i.e. for higher 
values of P, RV is smaller. The probability that the correlation coef- 
ficient is caused by random data is 7 per cent. Thus the correlation 
is not significant, however, one should keep in mind that we are 
clearly limited by low number statistics. As discussed above, this 
correlation is somewhat intrinsic as P is used for calculating RV. 

Further observations will have to show if the possible correla- 
tions between different parameters are spurious. 

Differences between the samples of FSRQs and BL Lac type 
objects are not significant due to the extremely small samples. Both 
the violence and the relative violence show similar ranges for the 
two samples. Different dependencies on redshift are not detectable 
as the redshift distributions of the two samples are dilferent. 



Figure 14. Dependency of relative violence RV on both redshift (Upper 
Panel) and degree of polarization (Lower Panel). Relative violence is in 
units of lO^^min"'. Non-variable polarized sources are also shown in the 
plot, for these we set the relative violence to RV = 0. However, we do not 
plot the values for the two possibly polarized RQQs. 



by the fact that we observe different restframe wavelength ranges 
due to the different redshifts of the objects. This would mean that 
variability is more violent for bluer restframe wavelength. How- 
ever, objects at higher redshift have to be intrinsically brighter so 
we might be observing a luminosity effect. 

To check for possible correlations between the violence V, rel- 
ative violence RV, degree of polarization P and redshift z we cal- 
culate the Spearman rank-order correlation coefficient for the data 
plotted in Fig.[T3]and Fig.[T4] The Spearman rank-order correlation 
coefficient is widely used in statistics to test how well an arbitrary 
monotone function describes the relation between two variables. 
We present the results of the Spearman test in Table|4] 

We calculate the Spearman-rank correlation coefficient p and 



6.2 Comparison of our results with those of other authors 

PINOV has not been studied systematically so far. However, several 
groups studied INOV and microvariability in RQQs and RLQs and 
tried to study if they show different behaviour concerning INOV. In 
this section we will thus sum up studies about microvariability and 
intranight variability in quasars so far and discuss the relations be- 
tween the observed INOV and PINOV properties of different types 
of AGN. 

iGopal -Krishna et alj ( Il993h studied two RQQs searching for 
micr ovariability but did no t detect any signs of the very. 

Ide Diego erZI ( ll998l) studied intranight variability in a well- 
selected sample of 34 quasars, half of them RQQs, half of them 
core-dominated RLQs. They found five out of 17 RQQ to show 
variability, for the RLQs, depending on the determination of vari- 
ability, 5-8 sources out of 17 showed variability. Using elaborate 



5 www.scipy.org 
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statistical tests. lde Diego et al. I( ll998l) did not find significant diff'er- 
ences in the observed variability between the two subsamples. This 
might i mply that INOV is not mainly dependent on radio-loudness. 

,Gopal-Krishna et al] ( 2000) observed 16 RQQs and one radio- 
weak quasar searching for microvariability. Determination of duty 
cycles is difficult for their sample as they used a wide variety of 
flags for the detection and non-detection of variability, owned to 
the fact t hat they tried t o specify the variability patterns observed. 

,Stalin et alJ dlOOSl) observed 20 high-luminosity QSOs, 12 of 
which were RQQs. INOV was detected in more than half of the 
objects. Collecting data from the literature they found that opposite 
to BL Lacs, RLQs and RQQs show similar and comparatively low 
duty cycles for INOV. This seems to imply that INOV is not mainly 
dependent on radio-loudness, but mostly on the orientation of the 
relat ivistic jet. 

ICarini et a i] ( l2007h observed a sample of seven quasars, five 
of them RQQs, two of them radio-intermediate quasars and did 
not detect intranight variability in any of them. They also collected 
data from the literature and found radio-loud objects to show much 
higher percentage of intranight variability than radio-quiet objects. 
However, contrary to other studies this sample is not well selected 
and different observing strategies make a clear decision about duty 
cycl es impossible . 

lGovaletai] ( l2007b monitored 1 1 RQQs over 19 nights for in- 
tranight variability and calculated a duty cycle of ~ 8-19 per cent 
using th e same statistic al limits for detection of variability advo- 
cated bv lde Diego et al.l |i998). 

A main problem comparing the different studies is the fact 
that different authors use different methods to decide if an ob- 
ject is variable or not. All kind of meth ods from ANOVA 
dde Diego et aljl993) over simple -methods l |jang|200^ to 'eye- 
ball decisions ' jGopal-Krishna et al.|[2000i) were used, producing a 
big amount of data tha t is diff icult to compare. 

lAndruchow et al.1 ( l2005h studied PINOV in a sample of 1 8 BL 
Lac type objects (eight of them X-ray selected, ten radio-selected). 
To our knowledge this is the only systematic study of PINOV in 
AGN published. They found high duty cycles for variability in both 
the degree of polarization and the position angle. Most of the ob- 
jects showed moderate to high degrees of polarization. However, 
their results can hardly be compared to our paper as a ^--criterion 
wa s used bv lAndruchow et"aL ( l2005h to decide about the variabil- 



ity. [p^^^^Rom^eUji 2006 ) presented a study in which system- 
atic errors for dual beam polarimeters were found. These system- 
atic errors are also known for ALFOSC0 and are immanent in dual 
beam polarimeters as used also bv Andruchow et al. (.2005). Thus 
we suspect that part of the detected variability in this study may be 
caused by these systematic errors . 

As lAndruchow et alj fcOOSi) , we find a extremely high duty 
cycle for polarized emission in BL Lacs. This is however not a 
surprise since polarized emission is one of the defining criteria of 
BL Lac objects. As for the duty cycles for PINO V in BL Lac type 
object s, we find similarly high duty cycles as in lAndruchow et all 
( I2OO5I) . Due to low number statistics the expected differences 
caused by the use of different statistical methods might not show. 

Studying not only microvariability but also polarization in this 
context clearly breaks the degeneracy between RQQs and RLQs. 
Their polarization properties are clearly distinct. This indicates that 
the emission mechanisms in RQQs and RLQs are clearly different. 



see |http://www.not.iac.es/instruments/alfosc/polarimetry/accuracy.html] 



6.3 Consistency with theoretical models 

Different theoretical models exist that explain variability on dif- 
ferent timescales. In this section we will present these models and 
discuss them with respect to our results. 

Czerny ( 2006) discussed the role of the accretion disk in AGN 
variability for a geometrically-thin, optically thick disk. The dy- 
namical time scale for such a disk can be calculated as follows: 



td\n 



(14) 



With a time scale of 4 h and a black hole mass of Mbh = 10' xM© 
this yields a radius of r < RschwarzscinU- Thus dynamical instabilities 
can not explain variability on intranight timescales. 
The thermal time scale can be calculated as: 



(15) 



with a being the disk parameter described bv lShakura & SunvaevI 
l ll973h . Assuming a characteristic value of a = 0.1, this again 
yields r < Rscinvarzschiu and thus thermal instabilities can not cause 
variability on intranight timescales. 

The viscous timescale can be calculated as: 



■ ^!henna! ^ 



(16) 



with hd being the disk thickness. So depending on the thickness of 
the disk, the viscous timescale could exhibit the right timescales to 
explain the intranight variability. 

Most theoretical models sadly lack timescale estima- 
tions, however, c ombinging theoretical models (see e.g. 
iMangalam & Wiit j [1993, Li & Cao 2008) with the time esti- 



mates from Czemv ( 2006h makes accretion disc instabilities quite 



unlikely on time scales of a few hours, only viscous instabilities 
can happen on reasonable radii for the given black hole mass 
and time scale. Thus it c an be stated that early assumptions by 
e.g. iGopal -Krishna et al.l ([l993) that accretion disk instabilities 
are observed as microvariability in RQQs, would point towards 
viscous disk instabilities. 

We interpret our results as follows: The dependency between 
high levels of polarization and PINOV seems to indicate that 
PINOV is a jet phenomenon as highly polarized emission is not 
expected from the accretion disk. This implies that physical con- 
ditions in the jet that cause PINOV have a very high duty cycle. 
Whatever causes these variations should be different from shock- 
in-jet phenom ena as described bvlMarscher & Ge^ j 19951) and ob- 
served e.g. bv lMarscher et alj 1 I2OO8I) as these shock fronts in jets 
evolve on timescales of weeks or months. However, turbulence in 
the shocks could explain PINOV and INOV. Changes in polariza- 
tion properties in blazars have been observed by e.g. Holmes et al] 
l ll984l) on timescales of several days. These kind of changes should 
be visible on intranight scale. However, physical conditions be- 
hind these observations are still unclear. Multiwavelength studies 
of such variations might be a way to distinguish between differ- 
ent models. Observations of this phenomenon so far have been 
described by several variable components in the j et with differ- 
ent magnetic field orienta tions and strengths (see e.g. lHolmes et al.l 
ll984lTakaloetalJll992h . However, this is not a physical descrip- 
tion. Theoretical models exist that can explain relatively fast vari- 
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ation in polarization p roperties (see e.g. Isiornsson & Blumenthall 
ll982llBi6ms"son|[ 1983) but have not been tested so far. 

Assuming jet emission as a main cause of PINOV implies that 
the jet emission in RQQs is either very weak or dominated by other 
regions of the jet than in radio-loud AGN. The low levels of polar- 
ization observed in this study might be due to residual background 
polarization or originate in scattering phenomena in the accretion 
disk or surrounding dust, however, this polarization does not seem 
to change on small time scales. The orientation of the jet also affects 
the brightness, variability and polarization properties. The smaller 
the viewing angle the more the radiation from shocks is beamed and 
thus the degree of polarization is higher, the objects seems brighter 
and variability appears to be faster. 

Sadly, we were not able to constrain flux intranight variability 
in our study, thus the link between flux and polarization intranight 
variability is still unclear. However, as radio-loud AGN show high 
duty cycles in both flux and polarization, they must coincide in 
many cases. Additionally, as we do not detect polarization variabil- 
ity in RQQs expect in one questionable case it can be assumed that 
RQQs show flux variability without varying in polarization. It is 
clear that we could explain pure intranight variability without po- 
larization variability by viscous instabilities in the accretion disk. 
We would not expect rapidly changing or high polarization from 
the accretion disk. Thus these kind of mechanisms are a promis- 
ing approach to explain intranight variability without accompany- 
ing polarization and PINOV. On the other hand, if we could affirm 
the connection between PINOV and INOV this would clearly point 
to jet phenomena causing the variability in radio-loud objects. 

In any case, there should be intermediate sources, thus study- 
ing PINOV and INOV in different classes of AGN should give a 
lead to the emission mechanism behind intranight variability. How- 
ever, extremely large samples are needed to overcome scatter in the 
properties of different samples. 



7 CONCLUSIONS 

We conducted one of the first systematic studies of polarization in- 
tranight variability. In this paper we present observations of a sam- 
ple of 28 quasars, 12 RQQs and 16 radio-loud blazars, at low to 
moderate redshifts (0.048 < z < 1.036) with a duration of 4 h each. 
We used statistical methods to detect polarization intranight vari- 
ability (PINOV) and measured the degree of polarization. 

We find that radio-loud and radio-quiet AGN show clearly dif- 
ferent behaviour. None of the RQQs show high levels of polar- 
ization. Only two out of 12 RQQs are polarized, both on levels 
P < \ per cent. On the other hand, 13 out of 16 radio-loud AGN 
show moderately to highly polarized emission. The three radio- 
loud AGN that are unpolarized are FSRQs, meaning that all BL 
Lac type AGN are polarized. Nine out of 15 radio-loud AGN show 
PINOV. For the FSRQs, four objects show PINOV, whereas one is 
polarized but non-variable. For the BL Lac type objects five show 
PINOV, three are polarized but non- variable. We find a potential de- 
pendency between degree of polarization and occurency of PINOV 
with all objects with P > 5 per cent showing PINOV, but not all 
variable objects showing degrees of polarization > 5 per cent. We 
do not find redshift dependencies of the measured degree of polar- 
ization, the occurrence of unpolarized emission or the occurrence 
of PINOV. However, our sample is relatively small, making the re- 
sults tentative. 

We detect two unpolarized objects showing polarization vari- 
ability, one of them is a RQQ, the other is a FSRQ. Such a be- 



haviour might either be explained by a flickering low level polar- 
ization or might be 'false'polarization variability. If there is under- 
lying constant but low polarized emission from e.g. scattering in 
the host galaxy or the dust torus, non-polarized varying flux from 
another component (e.g. the accretion disk) can cause variability in 

We study the violence of the variations in polarization. We do 
not find the violence of the variations to depend on the degree of 
polarization. We conclude that there might be a dividing line for 
variability at f > 5 per cent, with all objects showing high de- 
grees of polarization showing PINOV. However our sample is to 
small to detect a possible dependency between the violence and the 
degree of polarization. Another possible cause for the seemingly 
dividing line for PINOV might be the fact that detecting variability 
in low polarized objects is much more challenging. If we consider 
a highly polarized object with P = 25 per cent a comparatively 
small relative violence RV is easily detectable, whereas the same 
relative violence will not be detectable with the same error range at 
P X \ per cent. We detect a probable dependency between the vio- 
lence and the redshift. Still, this finding is not very definite due to 
small number statistics. A possible redshift dependency might be 
either a 'real'redshift dependency or e.g. wavelength dependency 
or might even be caused by selection effects as we observe intrinsi- 
cally brighter objects at higher redshift. 

From our results we conclude that PINOV and highly polar- 
ized emission originates from the jet. We conclude that the optical 
jet emission is either weak or emanating from different, less tur- 
bulent regions of the jet in RQQs. We also conclude that whatever 
causes polarized emission in RQQs is not able to produce high lev- 
els of polarization or rapid changes in polarization. The fact that 
PINOV is relatively common in blazars, especially when high de- 
grees of polarization are observed, poses the question which phys- 
ical conditions cause this phenomenon. In any case, mechanisms 
causing PINOV must have very high duty cycles in blazars. 

Future carefuU study of intranight variability in BL Lacs and 
FSRQs will have to answer the question what causes the variability. 
This will be a big step towards understanding the physical con- 
ditions in the jet. However, single filter observations might not 
be enough to reveal the physics behind PINOV. While variations 
on time scales of weeks and month are alre ady quite well under- 
stood and desc ribed in some cases (see e.g. iMarscher et alj|2003 . 
IValtonen et id]|2008i) , short-time polarization scales have been less 
studied so far. Another point that could not be addressed in this 
paper is the correlation between INOV and PINOV. However, we 
hope that this study emphasises the importance of understanding 
PINOV as it is very common in blazars. 
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